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I. INTRODUCTION 


Today, naval architects are greatly interested in the design of 
super tankers. However, many difficulties arise in design when the 
strength of the super tanker is considered from an economic point of 
view. The transverse strength problem is apparently one of these diffi- 
culties. 

There are several ways in analysing the transverse strength prob-~ 
lem. In each method, the two-dimensional approach is simpler, but the 
three-dimensional approach is so complicated that hardly no one wishes 
to use it. Therefore, some approximate is devised and sufficient results 
are achieved. But still some architects wish to know the three- 
dimensional effect on transverse strength analysis. 

Some prominent naval architects suggested solving the three- 
dimensional strength calculation but most of them did not apply numerical 
computations to existing complicated ships, for even though the computor 
is available, it is still no easy task. So many variables exist it is 
very tedious work to put them into computors. 

In this paper we are going to consider these tedious problems by 
using FORTRAN 7090. 

As previously mentioned, there are many ways to solve the three- 
dimensional transverse strength problem and each one has its own merits, 
although they will give almost the same level of accuracy. 

Since some approaches are so complicated and give the same level 


of accuracy as more simple approaches, they become useless. Therefore, 





we are going to deal with a well-known and more simple approach--the 
slope-deflection method--to attack these three-dimensional calculations 
and examine its effect. 

In large tankers especially, each span has such deep depth- 
breadth ratios that we suppose shear force contributes more than usual, 
bhepetore, We take into account shear force effect in slope deflection 
equation. 

Eventually, we will check our computed results by means of 
slope-deflection method with one of the fixed point methods which was 
given in our class by Professor Schade. This fixed point method differs 
in assumption from the slope-deflection method, but we can make this 
fixed point assumption comparable with our slope-deflection assumption 
by restricting our loading conditions. Although we are dealing with the 
fixed point method, it will not give us shear force effect and three- 
dimensional effect. We can see, therefore, only the differences between 
those two: the slope-deflection method with shear force effect and the 


fixed method without. They will lead us to some interesting conclusions. 


Il. BASIC SLOPE-DEFLECTION EQUATION CONSIDERED WITH SHEAR EFFECT 
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Using all of these results and the result of (1), we can easily get 


the following basic formula using the following notation. 
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the center of load for each span. 


For the purpose of our calculations, we only need two kinds of 


loading spans. 


They are shown on Table I and Table II. 
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Table 1 





Table 2 





III. SETTING UP SIMULTANEOUS EQUATION 

A. Assumption 

Once the designer has made the decision to consider the trans- 
verse frame of a ship, he is then faced with the task of deciding on 
simplifying assumptions. The main objective of these assumptions is to 
ereate a fictitious structure and condition of loading, which will re- 
duce the amount of labor in performing the structural analysis without 
compromising too severely the validity of the numerical answers. 

(1) One of the factors a designer must consider is the degree of 
dimensional similitude he can retain in the structure he selects for his 
computational work. For example, he must decide how much of the shell 
plating should be included in the analysis of a transverse frame as ef- 
fective material. We are going to use Effective Breadth Conception in 
Ship Structure Design, published in 1954 by Professor Schade. 

There are several other assumptions to be applied in this paper. 

(2) The relative angle initially formed by the tangents to the 
elastic curve at a point where two or more members are joined are 
assumed to remain unchanged after loading; that is, if two members are 
originally joined at right angles, the angle formed by the tangents to 
the elastic curve after the loads are applied will remain 90°. 

(3) Joints of 1-1 (or 9-2) and 4-2 (or 13-2) are fixed in space 
but free to rotate in three-dimensional analysis and additionally 9-1 
(or 14-2), 15-2 (or 13-1), 14-1 and 15-1 are fixed in space in two- 


dimensional analysis (see Fig. 2-2). 


B. Simultaneous Equation 
Before setting up simultaneous equations we must investigate the 


number of unknown variables. From the basic slope equation we can take 





up the following variables as unknown variables. 
1. 01(1) = 92(9) = X(1) 
2. @1(9 61(5) = e2(14) = x(2) 
3. @2(1) = 92(10) = 61(2) = X(3) 
nh. 92(5 61(10) = @1(6) = x(k) 


5. 92(2 62(11) = 61(3) = X(5) 
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13. ¢(3) = (7) = x(13) 
14. ¢(4) = 4(8) = x(14) 
15. 9(9) = g(10) = g(11) = g(12) = g(13) = x(15) 
16. ¢(14) = x(16) 
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Now we have 17 unknown variables in one transverse frame. 
Therefore we have to have at least Ll? simultaneous equations for each 
ice eg 

1. ML(1) + M2(9) = 


2. ML(9) + ML(5) + M2(14 
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3. M2(1) + M2(10) + M(2 


) 
) = 
4. M2(5) + ML(1O) + M1L(6) 
5. Me(2) + M2(11) + ML(3) = 
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6. M2(6) + ML(11) + ML(7 
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7. Me(3) + Me(12) + M(4) = 0 

8. M2(7) + mi(12) + M(8) = 0 

9. Me(4) + M2(13) = O 
10. Me2(8) + ML(13) + M2(15) = 
11. Q1(6) + Q1(2) - Q2(1) - Q2(5) = 0 
12. Ql(7) + Q1(3) - Q2(2) - Q2(6) = 0 
13. Q1(8) + QL(4) - Q2(7) - Q2(3) = 0 


14. X(11) x SPAN(1) + X(12) x SPAN(2) + X(13) x SPAN(3) + X(14) 
x SPAN(4) = 
For the three dimensional analysis we can set up the following 
three more equations. 
15. Yp= X(15) x SPAN(13) + X(17) x SPAN(15) 


Wes X(15) x SPAN(9) + X(16) x SPAN(14) 


an. 
| 


ive = X(15) x SPAN 
Here Lp » Ip and te are deflection of bottom and deck longitudi- 
nals and longitudinal BHD. We suppose that longitudinals are only acted 
by the concentrated load at joints with transverse frames. Then acting 


load Pp ; Py and P, are expressed as follows: 


PR = 2Q1(15) 
Pp = 2Q1(14) 
Pr = Q1(9) + Q1(10) + Q1(11) + Q1(12) + Q1(13) - Q2(15) - Qe(14) 


For the two dimensional analysis we need only to solve 14 simultaneous 
equations with unknown variables 14 by setting simply X(15), X(16) and 
X(17) to zero in the above equations. 

In this paper, for the practical numerical computation, we will 
take an existing Japanese built Tanker as an example. She has 47,000 
tons Dead Weight and one of the tank's length is 503.936"' and height 


604.329"* and its breadth is 1188.974"4 One of the web frame™® 
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scantling is shown in Fig. 1 and three web frames placed in a length of 


tank. 





FIG. 1 TRANS. WEB FRAME 





IV. PRELIMINARY CALCULATION 


A™™ Section modulus calculation 


From the effective breadth conception curve CASE (III), we can 


easily find out /#: Effectiveness ratio. From these data we can calcu- 


late moment of inertia and neutral axis. 


The results are as follows: 


MOMENT OF INERTIA 
OF THE SECTION INERTIA TO FREE FLANGE 
0.350 25839 .827 2.339 
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V. SOLVING A SIMULTANEOUS EQUATION USING COMPUTOR ‘/090 

Before solving a simultaneous equation we decide to solve two 
different loading conditions: (1) favorable loading condition (see Fig. 
3), (2) most critical loading condition (see Fig. 4). 

For the girder deflection we must also think about shear defor- 
mation because of having relatively short span length and deep depth. 
For this purpose we can find deflection equations of simple beam at the 
condition of both ends fixed and concentrated at some span point. 

One can use the next simple equation for the concentrated load 


at the center of the beam. 
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When the concentrated load is acting at the point of pee span 


length, deflection can be expressed as follows: 
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b>? P (1 + 2c_) 
p 512 p 


Ty 


Qo|\O 
ea] 


Where Bi Us res catieys at the point of load applied 
Ap.» Web area of b beam 
P ..-Concentrated lead 


For the computed deflection, see curves 4-1, 4-2 and 4-3. 
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FIG.4@1 UNIT LOAD DEFLECTION OF LONGITUDINAL CENTER GIRDER 
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FIG.4-2UNIT LOAD DEFLECTION OF DECK CENTER GIRDER 
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FIG. 43 UNIT LOAD DEFLECTION OF LONGITUDINAL BHD. 


10 


—_m ie owe ee 
—_ 


1/4 POINT LOAD 
—_ ———1/2 POINT LOAD 
_...--.....3/4 POINT LOAD 





Be 


One of the computed results is tabulated below as Table 5. 


Angle Deflection angle in Radian 
x aulors 


X(1) Cai fue 
Nae -0.413502 
X(3) -0.226263 
X(4) 0.128071 
X(5) 0.171848 
X(6) -0. 369586 
X(7) -0.501310 
X(8) 0.550122 
X(9) 0.647280 
ee) -0. 336503 
Kiel) 0.028208 
X(12) -0.080296 
X(13) -0.070333 
xe) 0.096957 


Table 5. Acquired unknown variables X at the 


forward frame (load 1) in two 


dimensional analysis. 


From these values we can get moment at edge 1 and 2 (left and 


right side or upper and lower side) for every span. 


the results in Table 6. 


You can see one of 
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SPAN No. MOMENT 1 x 102 MOMENT 2 x 103 
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Table 6. Moment of the forward frame 


From these end moments we can get very easily a middle portion 
moment by simply superimposing loading to each span. This work can be 
done very easily if we use the Computor for every span. The results 


are plotted in figures. 


VI. FIXED POINT METHOD 
A. General 
For the purpose of checking above results and investigating 


shear effect, we will apply here numerical computation by the method of 
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fixed point. The fixed point method is originated from the three moment 
equation. In this approach, we have the following assumptions: 

(1) All nodal points are fixed in space. 

(2) No shear force is considered. 
We will not discuss here in detail about this derivation (see Appendix), 
but will apply this method in practical calculation using the same load- 


ing conditions to the same scantling frame ship. 


Then the moment at the ends of the span can be expressed by the follow- 


ing formulas. 
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js is the area below the simple moment curve. 


and x is the distance from the left side to the center 


of the area Ag 


For our computation purpose we can find out the following data from the 


handbook. 
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handbook. 


hh, i" 


a 
© = 360. (8q,, + Ta, ) 


a 


seal (745 + Sq, ) 


cS 


I 


Re 





If once we get moment for each span, then we can compute each joint 
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moment using stiffness and fixed points by one-by-one adding process as 


follows: 





But in doing so, we must be very careful in the sign for some opening 


span. 


(Ref. illustrative calculation of moment distribution Table 12.) 


Now here is the practically computed results listed. 
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On the next page, f/f’ and A/A' tables are also prepared. 
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The computed end moment by its own span load are listed in 


Table 10. (We restricted our computation in the forward frame only.) 


LOAD I LOAD II 
SPAN No. 
MOMENT 1 x 10° | MOMENT 2 x 102 | MOMENT 1 x 102 | MOMENT 2 x 103 


0252657 ope atin Ze) 0.52657 





0.48144 C747568 O, tery 
On 7250 0.81965 0.79450 
1.03886 1.03658 1.03886 
-0.54524 
-0.49807 
-0.89650 
-0.13608 


-5 59423 -5 59423 


O O 


O -6. 82683 -6. 89544 


Table lO 


From the above computed end moment at edge 1 and 2 we can dis- 
tribute moment by using stiffness and fixed points P and A. (See 


Table 12.) Then we can get the following final end moment. 
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Table ll. Final end moment in tons-inch 


LOAD I eae 
SPAN No. 
MOMENT 1 x 102 | MOMENT 2 x 102] MOMENT 1 x 102 | MOMENT 2 x 103 
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Once we get each span's end moment then we can get various 
point moments very easily by superimposing process. And these numeri- 
cal computational results are plotted in figures. 

Also stress curves are available in this present paper, but we 
did not compute the moments and stresses for all the frames by the 
fixed point method. We only calculated the forward frame for both 
loading conditions. 

The reason is this: the fixed point method is, of course, a 
two dimensional one and each nodal point is fixed in space; therefore, 
a little difference of water head will bring us a very insignificant 


computational difference. 


VIL. DISCUSSION 

Most of the naval architects are anxious to learn about the 
three-dimensional effects on the transverse frame strength of the ship. 
But the computation is rather lengthy and tedious to work out so the 
architects estimate the effect on transverse frame strength by the 
Simplified method which satisfies them. In this paper we dealt with 
this lengthy calculation by use of the computor. It cannot be done 
without this machine because so many unknown variables are interacted. 

Even a casual glance at the graph will show that there are not 
many differences between the two-dimensional and the three-dimensional 
analysis. 

Especially in load condition @, the results are nearly the same. 
Also in load condition 2, we do not see much difference between these 
two. 

In load condition 1, at the forward frame (Fig. 5) we can see 


also almost no differences between the fixed point method and our 
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present method in the upper half of the frame that is span No. l, 2, 5, 
6, 9, 10, 11 and 14. That means that lightly loaded spans do not have 
any shear effect at all, but loading is steadily increased, then the 
discrepancy is noticeable. 

We are now going to examine these graphs individually for more 
detail. 

A. Load Condition 1 

(1) This@lipadine condaitien 1s ballasted withesea water in 
the center tank and the side shell is exhibited to the hydrostatic load. 
This is one of the most favorable loading conditions. In this loading 
condition, the forward frame does not show us much discrepancy between 
two-dimensional and three-dimensional analysis. In this loading condi- 
tion, junction points 1-2 and 2-1; 2-2 and 3-1; 3-2 and 4-1; also 5-2 
and 6-1; 6-2 and 7-1; and 7-2 and 8-1 remain in space because of the 
Symetric loading condition. Therefore the fixed point method and our 
Slope-deflection method can be comparable in quantities. 

Thus, this means if we neglect shear effect, we will lose al- 
most 4,000 Tons-inch moment at the span edge which is most heavily 
loaded. 

(2) Middle frame: We can see a little difference between 
these two methods, but not too much. This means that the middle frame 
is the most effected frame among the three frames because the girders 
will deflect most at the center. 

(3) In the aft frame, we can see almost no effect of the 
three-dimensional effect on the moment curve. Almost the same as the 


forward frame. 
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Be Load Condi tiom 2 

(1) We think this loading condition is most critical, but 
still two-dimensional analysis and three-dimensional analysis have the 
same aspect. They don't change very much, but here are several inter- 
esting aspects being displayed. That is, in our method we fixed joint 
1-1 and 9-2 and 4-2 and 13-2 so that as a whole side frame and BHD bent 
something like a semi-fixed beam at the bottom and deck, but by the 
fixed point method nodal points are fixed in space so that they have 
almost the same aspect as in load condition 1. But from this curve we 
can see also a big difference in heavily loaded lower spans. 

(2) In the middle frame we can see some conceivable dis- 
crepancies between those two methods, but still in good shape; in 
three-dimensional analysis, however, the moment at point 15-1 and 9-2 
is reversed in direction because of the restraints by girders released 
at 15-1 and 13-1 points. The other spans have almost the same results, 
nevertheless one can see a little higher moment in the frames 9, 10, ll, 
i2eancde 1. 

(3) In the aft frame, the aspect is almost the same as in 
the feneda frame. We cannot see any noticeable discrepancies between 


those two approaches. 


CONCLUSION 

Most designers of super tankers are worried about the three- 
dimensional effect on the transverse frame strength. But from our com- 
putation we cannot find large discrepancies between two-dimensional and 
three-dimensional analysis. However, in the three-dimensional analysis, 


calculation is more detailed and tedious than in the two-dimensional. 
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For our present discussion, three-dimensional computation is 
more academic than practical, otherwise we would not use any simplified 
design process such as Professor Schade's Design Curves for Cross 
otiffened Plating. 

Thus our conclusion is: forget about the three-dimensional 
analysis, unless you intend to use any simplified design process pre- 
pared for the three-dimensional analysis; but take into account the 


shear effect on the two-dimensional analysis on super tanker designs. 
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FIG. 5 FORWARD FRAME MOMENT CURVE (LOAD 1) 
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FIG. 6 MIDDLE FRAME MOMENT CUPVE (LOAD 1) 


TWO DIMENSIONAL 
__. __ _. THREE DIMENS TIONAL 


\ 


el 
——_a 





___2000 IN-TONS 








Sek: 


FIG. 7 AFT FRAME MOMENT CURVE (LOAD 1) 
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FIG. 8 FORWARD FRAME MOMENT CURVE CLOAD 2) 


TWO DIMENSIONAL 
— —_ — -[A Ree DIMENS TONAL 
See ee eo LNT ht T HOD 








ho. 


FIG. 9 MIDDLE FRAME MOMENT CURVE (LOAD 2) 
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FIG.10 AFT FRAME MOMENT CURVE (LOAD 2) 
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FIG. LL FORWARD FRAME STRESS CURVE (LOAD 1) 
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FIG. 12 MIDDLE FRAME STRESS CURVE (LOAD 1) 
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FIG. 13 AFT FRAME STRESS CURVE (LOAD 1) 
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FIG. 14 FORWARD FRAME STRESS CURVE (LOAD 2) 
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FIG. 15 MIDDLE FRAME STRESS CURVE (LOAD 2) 
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FIG. 16 AFT FRAME STRESS CURVE (LOAD 2) 
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IX. APPENDIX 


1. The three moment equation 
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Then by using the second moment area theorem to obtain 
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simplifying the above expression, the three moment equation is: 
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2. Fixed point method 





By three moment theory (Ref. Appendix 1) first three support (from left) 
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Apply three moment equation to span C and D 
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In general 
An 
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OTE: The position Of zero moment 1S 92 Tunction of geometry and the 
degree of fixity at the ends and independent of the type of loads or 


the magnitude of loads. 


LOADED SPAN E 





Apply three moment equation to span D and E (D is unloaded, E is 


loaded). 
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Where R, = Moment of field load moment area about right 
Support of span EL 


L. = Moment of field load moment area about left 


support of span E 
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stiffness factor is defined as the reciprocal of the angle of 
rotation caused by a unit moment. 
From the elementary mechanics of materials we know the 


following relation between angle of rotation and moments 
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By these relations we can easily obtain 
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If we use the notation ¢& _ to denote the stiffness of the left 


hand side member, and substitute ve Dye Ce ) we will get 
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In general, the stiffness of the left hand member at the right 
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In general cl 
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If we use the notation %& including left and right hand 


members! stiffness, then we have 
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oimilarly we can apply the above formula from right side to 


left side. The results are: 
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